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Abstract: We investigate spin correlation effects in the "antler" event topology pp —> A — )• 
B\Bi — > (l~Ci){l + C2) at the LHC. We study the shapes of several kinematic variables, 
including the relative pseudorapidity, relative azimuthal angle and the energies of the two 
leptons, as well as several mass variables M#, M e g, v^Wn' Mt2, Mct and Mqtx- We fo- 
cus on the two kinematic extremes of yfs — threshold and infinity — and derive analytical 
expressions for the differential distributions of several variables, most notably the cos#|_^ + 
variable proposed by Barr in hep-ph/0511115. For all possible spin assignments of particles 
A, B and C, we derive the eos#| _« + differential distribution at threshold, including the ef- 
fects of spin correlations. Our analytical results help identify the problematic cases for spin 
discrimination. 
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1. Introduction and motivation 

Once a signal of new physics is discovered at the Large Hadron Collider (LHC) at CERN, it 
must be interpreted in terms of the production of new particles with definite properties: mass, 
spin, charge, coupling strength, chirality, etc. Typically, one studies the properties of new, 
short-lived resonances by examining their decay products which are seen in the detector. For 
example, the electric charge of the parent particle (a short-lived resonance) is easily obtained 
by adding up the electric charges of its daughters (particles from a decaying process) , the mass 
is similarly determined by forming the combined invariant mass of the daughter particles, etc. 

Notoriously, this procedure breaks down when some of the daughter particles are very 
weakly interacting particles and are not visible in the detector (the prototypical examples 
being the neutrinos of the Standard Model (SM) or hypothetical dark matter WIMPs). Mass 
and spin measurements then become very challenging, especially at hadron colliders where 
the partonic center-of-mass (CM) energy and total longitudinal momentum in the event are 
a priori unknown. This has prompted a long line of research into developing new kinematic 
methods for determining the masses (see, e.g. [1,2] and references therein) and spins [3-34] 
of new semi-invisibly decaying parent particles. 

In this paper we tackle the more challenging of those two issues, namely, the determina- 
tion of the intrinsic spin of the new particles. Broadly speaking, most methods proposed so 
far in the literature fall into the following categories: 

• Methods applicable to the case of lepton colliders, where the CM energy and longitu- 
dinal momentum of the initial state are known [4, 17]. Here we shall be interested in 
hadron colliders like the LHC, where those methods would not be directly applicable. 

• Methods attempting either exact [27] or approximate [21] reconstruction of the kine- 
matics of the invisible particles on an event by event basis. Unfortunately, in order to 
be able to solve for the missing momenta, one needs a relatively long decay chain, with 
a sufficiently large number of intermediate resonances [35]. In contrast, in this paper we 
shall consider the case of the shortest decay chain possible (a chain with a single two- 
body decay), for which there are not enough kinematic constraints for reconstructing 
the missing momenta, so that one has to resort to likelihood methods [26,29]. 
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• Methods which study the shapes of the invariant mass distributions of the visible parti- 
cles originating from one given decay chain [3,5,6,11,12]. These methods again require 
a relatively long decay chain, which would produce at least two visible SM particles, 
and are not applicable to the short decay chain which will be considered here. 

• Methods which simultaneously utilize the measured visible particles from both decay 
chains 1 . Here one studies the distributions of suitably defined quantities, e.g. angular 
variables [8], asymmetry-like variables [29,33], or global event variables like M c g or 
fhr [30]. 

In this paper we shall focus on the most challenging case for spin determination: that of 
a single stage, two-body decay chain, in which a new parent particle B is produced and then 
decays to a single visible SM particle c and an invisible daughter particle C: 

B^c + C. (1.1) 

Motivated by the dark matter problem, we shall further assume that the new particles C and 
B are Z 2 -parity odd, so that the LHC produces the parent particles in pairs: 

pp ->■ BB — >■ ccCC. (1.2) 

In general, the process ( |1.2|) can be due to s, t and u channel diagrams. Which of those 
diagrams are present and/or dominant is a very model-dependent issue, which also depends 
on the nature of the visible particle c. For example, if c carries color (i.e. it is a quark or a 
gluon), one generically expects all three types of diagrams to be present. On the other hand, 
if c is a SM lepton, which carries lepton number but no color, then it is reasonable to expect 
( |l.2j ) to be mediated by s-channel diagrams only 2 . In what follows, for simplicity we shall 
assume that the SM particle c is a lepton £, and that the process ( |1.2|) is due to the exchange 
of a single ^-parity even particle A in an s-channel. The resulting event topology is shown 
in Fig. m and has been dubbed "the antler event topology" [36], due to its resemblance to the 
headgear of the common wildlife species in Wisconsin. In general, the s-channel particle A in 
Fig. D can be a SM gauge boson (a photon or a Z-boson), a SM Higgs boson, or some other 
new bosonic particle. Depending on the relative masses of particles A and B, the s-channel 
resonance A could be on-shell (if Ma > 2Mb) or off-shell (if Ma < 2Mb)- In the following 
we shall discuss both of those cases. 

In summary, as illustrated in Fig. jj], the process under study here will be 

pp A -> B 1 B 2 -»■ £ + rCiC 2 . (1.3) 

1 In a typical new physics model, the longevity of the dark matter particle is enforced by an exact sym- 
metry, the simplest option being Z2 parity, under which all SM particles are even, while the new particles 
(superpartners, KK-partners, etc.) are odd. A consequence of such Z2 parity is that the new particles are 
pair-produced, and each event has two cascade decays. 

2 In principle, t and u channel diagrams cannot be ruled out completely — for example, the dark matter 
particle C itself may also carry lepton number (although the sneutrino option in supersymmetry is disfavored 
by experiment), or the t and u channel diagrams could be mediated by exotic particles carrying both color 
and lepton number (e.g. leptoquarks), or the process (|l.2j) may involve lepton number violating interactions. 
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Figure 1: Left: the antler event topology under consideration in this paper. Bi and Cj are new i?2-odd 
particles, whose spins (as well as the chirality /? and 7 of their couplings) are a priori unknown. The 
intermediate resonance A is a i?2-even particle which may or may not be on-shell, and furthermore, its 
spin (and the chirality a of its coupling to the initial state partons) is also unknown. Right: pictorial 
definition of the kinematic variables used in this paper. Starred quantities refer to the CM frame of 
the B1B2 system while primed variables refer to the individual rest frames of the parent particles Bi. 

An important benefit of considering the leptonic (i.e. c = £) version of the antler topology is 
that it does not suffer from large QCD backgrounds (as would be the case if c were a jet). 
Another advantage is that it allows us to safely focus on the s-channel topology only, where 
we only need to introduce the very minimum number of unknown parameters (see Section |2] 
below) and tackle the problem of spin measurements in a rather model independent way. 

In this paper, we shall analyze the kinematics of the antler event topology (|T^) with 
several goals in mind: 

• First, we shall perform a comparative study of the several kinematic variables already 
discussed in the literature, supplementing them with a few variables of our own (see 
Sec. H). We shall be mostly interested in the question to what extent the shapes of the 
different kinematic distributions are sensitive to spin effects. The main objective of this 
initial study (presented below in Section |j) will be to identify the kinematic variables 
which retain the largest amount of spin information, and are thus most suitable for 
measuring the spins of particles A, B and C. 

• Next, we shall address the question how to infer the presence of spin effects in the distri- 
bution of any given kinematic variable. Clearly, in order to firmly establish the presence 
of spin effects, one must first know the expected shape of the distribution in the absence 
of any spin correlations, i.e. in the pure "phase space approximation". Then, any devi- 
ation from the expected "phase space" distribution will indicate spin effects, although 
a priori it will not be immediately clear whether those are due to non-trivial intrinsic 
spins of the particles {A, B,C}, non-trivial coupling chiralities, or some combination 



- 4 - 



of the two [20] . In Section |B| we therefore provide the expected "phase space" distribu- 
tions for several kinematic variables, including the angular variable cos0^_^+ proposed 



in [8] (see Sec. 3A below for its exact definition). We shall pay special attention to the 
dependence of the kinematic variables on the total CM energy yfs in the event, and in 
Section |5,1| we shall derive a new formula for the shape of the cos «+ distribution at 
threshold (y/a = 2M B ). 

• Having observed a deviation from the expected "pure phase space" kinematic shape, 
one can be sure that non-trivial spin effects are in play. But how much of the effect is 
due to the spins and how much is due to the chirality structure of the couplings? This 
question is tackled next in Section ||. In order to explicitly separate the two effects, we 
derive analytically the cos0L» + distributions at threshold for various spin assignments 
and arbitrary chiralities. These results provide useful intuition and guidance for the 
study which comes next — the measurement of the spins of particles A, B and C. 

• Section ^ contrasts the shapes of the cos 0^_ £+ distributions obtained for various spin and 
chirality assignments. The ultimate goal there is to identify which spin combinations 
for the particles {^4, B,C} can be distinguished (at least in principle) in such a model- 
independent way (i.e. without any additional assumptions about the chirality of the 
fermion couplings) and which can be confused with each other. We shall specify the 
dangerous "twin" spin scenarios with similar cos#^_^ + distributions and discuss the 
possibility to discriminate them through the distributions of other, complementary 
kinematic variables. 

Before taking up these four goals, we first introduce our setup and notation in Section |2[ 
2. Setup, notation and conventions 

In this paper we study the process ( |1.3| ) depicted in Fig. |l[ The two particles C\ and C2 are 
invisible in the detector, while the two massless leptons £ + and t~ are visible. The double 
lines in the antler event diagram of Fig. |l| denote particles whose spin is unknown and needs 
to be determined. The masses of the particles A, Bi and Cj are also initially unknown, 
but can eventually be determined from purely kinematic methods [35-47], in spite of the 
challenge presented by the short decay chain. In what follows, we shall therefore assume 
that the masses Ma, Mb and Mq of the new particles have already been measured and we 
concentrate on the issue of spin determination alone. 

In order to perform a model independent spin measurement, we will need to separate 



pure kinematic effects (see Section [2.1| ) from spin and chirality variables (see Sections 2.2 and 



2.1 Kinematics 

In what follows, we shall find it convenient to describe the kinematics of the process ( |1.3[) in 
terms of variables which refer to three different reference frames (two of those are illustrated 
in Fig. Ill): 
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1. The LAB frame. This is the frame of the experiment, in which the visible particles 
are measured. Variables related to the LAB frame will carry no special designations, 
e.g. the 4-momenta of the two visible particles (the leptons) will simply be 



P*_ = Pf = [E 1: P 1T ,P lz ) =(E 1 ,P 1 ,Q 1 ,i Pl ); (2.1) 
P£ = p% = (e 2 , P 2T , P 2z ) = (E 2 , P 2 , 2 , <p 2 ) , (2.2) 



where 0j and <pi are the usual spherical angular coordinates. Note that the angles 
Oi are measured from the beam axis (z-axis) and ipi are measured in the fixed plane 
(common to all events) transverse to the beam. Since the leptons are assumed to be 
massless, we also have Pi = E^. The LAB variables are not shown in Fig. |]. 

CMBB: the CM frame of the parent pair B\B 2 . This is the relevant frame for 
describing the 2 — > 2 process pp —> B\B 2 and the corresponding kinematic variables 
will be denoted with an asterisk (Fig. [l]). The final state of the B\B 2 pair is described 
by three degrees of freedom, which can be taken to be: the CM energy ys* and the 
spherical angular coordinates of particle B\, which are 0* (still measured from the 
beam axis) and tp*. Since y/s* and <p* are invariant under longitudinal boosts, they are 
the same in the LAB frame: 



s* = yfs, (2.3) 
<P* = <P, (2-4) 

so that in what follows for simplicity we shall omit the asterisks on y/s and (p. Instead 
of y/s, one could equivalently consider the (magnitude of the) momentum P* of the Bi 
parent particle 

or the corresponding boost factor rf 

rf = cosh" 1 (-^T-l = sinh_1 • ( 2 - 6 ) 

' \2M B ) \M B ) y ' 

The main variables P* , Q* and <p>* of the CMBB frame are depicted in Fig. |l| in blue. 

3. CMB1 (CMB2): the CM frame of an individual parent Bi. The kinematics of 
each final state lepton is most easily described in the rest frame of its parent particle, 
where variables will be denoted by a "prime" superscript (Fig. |l|) and will also carry 
a subscript "1" or "2" to indicate the CM frame of B x CMB1 or the CM frame of B 2 
CMB2. The three degrees of freedom describing the kinematics of each lepton can be 
taken to be: the magnitude of the momentum P?, the polar angle 9^ and the azimuthal 
angle Note that, as shown in Fig. |l|, the polar angles 9[ are measured with respect 
to the direction of the parent particle Pi and not the beam axis, as was the case for the 
angle 0. Similarly, the azimuthal angles are measured differently from the azimuthal 
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angle ip in the CMBB frame. The momentum of each lepton in the corresponding parent 
CM frame (CMB1 or CMB2) is a constant which only depends on the mass spectrum: 



In place of the polar angle 0\ one could alternatively use the corresponding pseudora- 
pidity 



Vi = - In 



tanh 



(2.8 



The frame CMB1 (CMB2) is obtained from the frame CMBB by boosting with a Lorentz 
factor 7 = cosh rf along the direction of B\ (B2). 

2.2 Spin assignments 

Since the particles A, B and C are not directly measured, their spins are a priori unknown, 
and each can be a spin zero scalar (S), a spin ^ fermion (F), a spin 1 vector particle (V), etc. 
Following [11,20], we shall consider all allowed spin assignments 

(ABC) € j(SSF), (SVF), (VSF), (VVF), (SFS), (VFS), (SFV), (VFV) j (2.9) 

in the antler event topology of Fig. |l] with spins up to one, i.e. consistent with renormalizable 
interactions. The Feynman diagrams for these 8 cases are listed in Table ffl, together with a 
representative example from a popular new physics model. 

For the underlying 2—7-2 process pp — > BB, there are in principle 6 possibilities for the 
spins of the (AB) pair: 

(AB) e ((SS),(SV),(VS),(W),(SF),(VF)1 (2.10) 



which will be discussed in Sections 3.2 and 6.3. 



2.3 Coupling assignments 

As we shall see below, the shapes of the kinematic distributions are affected not just by the 
spins of the intermediate particles, but also by the chiralities of their couplings. Since the 
observed particles (leptons) are fermions, the corresponding vertices of the antler diagram may 
contain chiral projectors. In the spirit of Refs. [20,34], in each vertex of the antler diagram, 
we allow for an arbitrary linear combination of right-handed and left-handed couplings. For 
example, when C is a fermion and B is a vector boson, the BtC coupling is given by 

Hint = ^ily. (clPl + c R P R ) 4> c A b , (2.11) 

with arbitrary coefficients cl and cr. As indicated in Fig. [l], the relative chirality of that 
vertex will then be parameterized in terms of a single parameter 



r 2 



4 + 4 



(2.12) 



-7- 



Spin assignments 


Antler diagram 


Examples 


SSF 


> 6 




SVF 


> 4 


H° W + W £+£~vv 


SFS 




H 2 -> l+q -> £+£-B H1 B H1 [48] 


SFV 


y< 

/ \ 


H 2 -»■ -> £ + £-B 1 B 1 [48] 


VSF 


/ \ 


Z' -> £+r -> ^+^-x?Xi [36,49] 


VVF 


w 

/ \ 


Z' VF+TF" f+^P [50] 


VFS 




Z 2 -> -> £ + £-B H1 B m [51,52] 


VFV 


H 


Z 2 -)• -»• £ + £-B 1 B 1 [6,53] 



Table 1: The 8 different spin assignments for the unknown particles A, B and C in the antler 
topology of Fig. [l| The last column lists a few typical examples from popular theoretical models, 
where supersymmetric particles are denoted by a tilde, while a subscript n indicates a level n Kaluza- 
Klein mode in theories with Universal Extra Dimensions. 



so that 7 = corresponds to a vector-like coupling, while 7 = ±1 is a purely chiral coupling. 

One can extend this technology to the other two vertices as well [20] and define the 
relative chirality of the AB\B 2 vertex as 



and the relative chirality of the qqA vertex as 



a 2 
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Pos. Vertex 



Lagrangian term 



Vertex 



a: / (uP q (P L + e i5 «P R )i{j q a(P L + e iS «P R ) 

~ \ 

/ ijjqj/j, (a L P L + a R P R ) ip q A^ 7^ (a L P L + a R P R ) 

P: \ ^B^(bLP L + b R P R )il)BA^ ln(b L PL + b R P R ) 

b TpB (P L + e i5 »P R ) ^ B ct>A b [P L + e iS »P R ) 

\ F^FV b ((p - kYgT + (q- p) u g pp + (k - qYg pv ) 



(D^ B y(D^ B ) b(p-k)P 
bA B A^ B <t> A bg^ 



bcf> B (j) B (j)A b 




7 : i ^Hii (clPl + c R P R ) ^ B A P C + h.c. 7 M {c L P L + c R P R ) 

-J 

\ 

i\>i (c l Pl + c R P R ) ipB<f>c + h.c. (clPl + c R P R ) 

\ ipllij, (°lPl + c R P R ) ipcA B + h.c. 7 M (c L P L + crPr) 



ipl (clPl + crPr) ipc4>B + h.c. (c l Pl + crPr) 



Table 2: Summary of the Lagrangian terms and the corresponding Feynman rules for the vertices in 
the antler diagram. 

Our notation is chosen so that the chiral coefficients ai, a R and the relative chirality a refer 
to the production vertex qqA of particle A; b^, b R and the relative chirality 3 j3 refer to the 
production vertex AB\B<i of particle B, while cl, c r and 7 refer to the production vertex 
BIC of particle C. 

Table lists the explicit form of the coupling in each vertex for all possible spin assign- 



? Of course, the definition of /3 does not apply when A and B are both bosons. 
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ments for A, B and C which we are considering. 



3. Kinematic variables for the antler topology 



Particle 


SUSY-like 


UED-like 


A 


7*/Z* 




7*/Z* 




B 


Zr 


176.62 


h 


505.45 


C 


X? 


158.18 


7i 


500.89 



Table 3: Mass spectrum in GeV for the two study 
points considered in the text. The SUSY-like mass 
spectrum corresponds to the CMS LM6 study point 
[54], while the UED-like spectrum comes from the 
Minimal UED model with i?" 1 = 500 GeV, XR = 20 
and m h = 120 GeV [4-6]. 



At a hadron collider like the LHC, the 
energy and momentum of the initial state 
partons vary from one event to another. 
The lack of knowledge about the kinemat- 
ics of the initial state prevents us from com- 
puting the momenta of the missing parti- 
cles Ci in the simple antler topology. (At 
a linear collider, where the energy and mo- 
mentum of the initial state are known, the 
missing momenta in the antler topology can 
be reconstructed only up to a two-fold am- 
biguity [17], even assuming precise knowl- 
edge of the masses of the intermediate par- 



ticles). Thus we are limited to variables constructed out of the measured 4-momenta 
(defined in eqs. ( |2.1| - |2.2D ) of the visible particles in the LAB frame. In this section, we 
review a number of such variables which have been discussed in the literature. In order to 
gauge their sensitivity to spin effects in the antler topology, we shall investigate their distri- 
butions for two conventional study points, whose mass spectra are listed in Table ||. For each 
of these two mass spectra, we shall consider two spin configurations: VSF for SUSY and VFV 
for UED, postponing the comprehensive analysis including the remaining spin combinations 
until Sections ^ and [7|. The resulting four cases will be contrasted in Fig. pi where we shall 
plot the unit normalized distributions for all variables introduced in this section. 



3.1 The relative pseudorapidity variable cos#J_» + 

In [8] A. Barr introduced an angular variable 0^_ f+ which is related to the relative pseudora- 
pidity of the visible particles in the LAB frame as 



cos ( 



1-1+ 



tanh 



m- 



Ar]f,- i+ 



%+ =vi- m- 



(3-1) 
(3-2) 



For brevity, from now on we shall use Cb to denote cos • The pseudorapidities of the 
two visible particles in the LAB frame are defined as 



1 



In 



Pi + Pi 



2 \ Pi Pi z , 

in terms of the measured momenta Pj and Pi z from eqs. (2.1-2^2]) 



(3-3) 
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3.2 Azimuthal angular difference /S.ip 

The azimuthal angular difference of the visible particles 



has also been shown to be sensitive to spin effects [28, 29] and we shall include it among our 
set of observables here as well. 

3.3 Invariant mass Mg+g- 

One can also consider the invariant mass of the two visible particles in the antler topology 



M i+l - = J(P X + P 2 r (Pi + P 2 ) u . (3.5) 



i-i 

The shape of the Mg+g- distribution was studied in detail in [36] for the case of resonant 
production (Ma > 2Mb)- It was shown that the Mg+g~ distribution exhibits a point (called 
a "cusp"), where it is not continuously differentiable. If the mass spectrum is such that 

M A < (2Mb) cosh ra (2.26) Mb, (3.6) 

then the cusp appears at the very peak of the distribution. On the other hand, when Ma > 
2.26Mb, the peak is smooth and the cusp appears somewhere below the peak [36]. 

In what follows, we shall prefer to work with dimensionless variables, so we need to 
rescale our mass variables by a suitable mass constant. Since the mass spectrum of the antler 
topology is assumed to be already known, the momentum Pi of each lepton in its mother's 



rest frame (2.7) is also known and we can use it for normalization. We therefore define the 



rescaled invariant mass 4 as 

M e+ g- = • (3-7) 

3.4 Contransverse mass Mqt 

The contransverse mass variable 



M CT = \2[ P 1T P 2T + P 1T ■ P 2T (3.8) 



proposed in [55] is also suitable for the antler topology. It is invariant under longitudinal 
boosts from the LAB frame to the CMBB frame, and under "back to back" boosts of the 
parent particles. The Mqt distribution has an upper kinematic endpoint 

M^ = 2P{ = M B (l - ^f) > (3-9) 



4 Rescaled quantities will be denoted with a hat. 
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which can be used for mass measurements [55]. Here we shall focus on the shape of the 
distribution of the corresponding rescaled Mqt variable 



Mor = ^. (3.10) 



Eqs. (|3.9| ) and ( 3.1C| ) imply that the endpoint of the Mct distribution is always at Mct = 1 



and this will be evident in the plots below. 

3.5 Doubly projected contransverse mass Mqtx 

The doubly projected contransverse mass variable Mct± was proposed in order to remove the 
upstream boost effect from initial state radiation (ISR) [41]. Mct±_ is the one-dimensional 
analogue of (|3.8| ) which uses the transverse momenta components which are orthogonal (hence 
the subscript "_L") to the upstream which can be caused by ISR or previous decays. Here 
we shall not consider ISR effects, so we are free to choose any one of the two transverse axes 
to define the "_L" components, e.g. 



M CTx = V2 (\Pi x P 2x \ + Pi A) • (3.11) 
The corresponding rescaled variable is 

Mct. = (3-12) 

and its endpoint is also at Mqtx = 1- F° r our purposes, Mqtx has an important advantage 
— the analytical formula for its differential distribution (in the absence of spin correlations) is 
already known [41], providing an important benchmark for spin studies (see Sec. |Oj below). 

3.6 Effective mass M e g 

The effective mass variable M e fj [56] has been widely utilized in SUSY searches. When applied 
to the antler topology, it reads 

M eff = Pit + P 2 t+ ■ (3.13) 

Here we use the rescaled variable 

M eff = ^- (3-14) 

Many SUSY analyses use M e ff as one of their selection cuts, since it is correlated with the 
mass scale of the particles produced in the event. 

3.7 The \/s mm variable 

The \/s miri variable was advertized in [57, 58] as a better estimator of the mass scale of the 
hard scattering, since it explicitly accounts for the masses of the invisible particles, and is 
defined in a theoretically rigorous way. y^mm ^ s the minimum value of y/s which is required 
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in order to account for the observed visible particles and the measured Et in the event. In 
the case of the antler event topology, the general formula for v^mm fr° m [57] reduces to 



V~S min = V M iV + \ P IT + P2T\ 2 + V 4M C+ B T > ( 3 - 15 ) 

already accounting for the fact that the combined mass of all invisible particles present in the 
event (in this case the two C's) is 2Mc- It was noted in [57] that for typical SUSY events the 
\/s min distribution peaks near the mass threshold (2Mb) for production of the two parents. 
This is why we choose to rescale the \/s min variable by 2Mb and define 

■ r - - ^ min (3.16) 



2Mb 



which we expect to peak near v^min — 1? as confirmed in Fig. 0. 

We anticipate that the y/s min variable will be useful for spin studies, particularly in the 
off-shell scenario when Ma < 2Mb- In that case, the two B particles are produced near 
threshold, and the threshold suppression of the yfs distribution will be sensitive to the spins 
of the A and B particles (see Sec. [T^ below). Since y/s m - m is designed to correlate with the 
actual yfs in the event, any spin effects present in the ^fs distribution will be inherited to 
some extent by the ^ min distribution as well. 

3.8 Lepton energy Ef 

The lepton energy Ei is a classic variable for kinematic studies at a linear collider. Here we 
revisit this variable for the case of a hadron collider, rescaling as usual 

Ei = § ■ (3.17) 

In what follows we shall not distinguish the two lepton energies E\ and E% and will always 
show them together on the same plot as E(_. 

3.9 The Cambridge Mt2 variable 

The Cambridge Mt2 variable was originally proposed [59] for precisely the antler topology 
case of Fig. Q. Since then, it has been extensively applied for measuring the mass spectrum 
of the new particles A, B and C. In the case of the antler topology, Mt2 is given by [42,60] 



M T2 (M C ) = v At + y At + M£, (3.18) 
A T = - (Pi T P2T + Pit ■ P 2 t) . (3.19) 



Just like \/s mm , Mt2 requires an input mass for the missing particle, which in eq. (|3.18|) has 
been taken to be the correct value Mq- In that case, the Mt2 endpoint is equal to the parent 
mass Mb- 

M^ ax) (M c ) = max{M T2 (M c )} = M B - (3.20) 



- 13 - 



This motivates us to rescale the M?2 variable as 

jCr _ Mt2(m c )-m c 

Mt2 ~ M B -M C ' (3 - 21) 
so that the distribution of Mt2 ranges from to 1. 

4. Comparison of different variables 

Having defined the nine kinematic variables of interest in the previous section, we shall 
now compare their sensitivity to spin effects in the antler topology. In order to get some 
preliminary idea about this, in Fig. ^ we show the distributions of our nine variables in the 
four cases, obtained by pairing the two mass spectra from Table |3| (SUSY-like or UED-like) 
with the corresponding spin assignments (VSF for the MSSM and VFV for MUED). If a 
particular distribution is affected by spins, one would expect a visible difference between 
the two dotted lines, showing the two different spin scenarios for the same SUSY-like mass 
spectrum. Similarly, one would also expect a difference between the two solid lines, which 
show the two different spin scenarios for a UED-like mass spectrum. 

Fig. [2] reveals that the nine variables from Section || exhibit varying degrees of sensitivity 
to spins. The distributions of some variables, like Mqt, Mctx and M?2, show very little 
variation and appear to be relatively insensitive to changes in the spins and masses. The 
case of Mqtx is particularly striking — all four Mqtx distributions are virtually identical. 
This begs the question whether the Mqtx distribution has any dependence on the spins at 



all. This issue will be tackled in Section 5.4 below. 

At the same time, Fig. |2| also demonstrates that the remaining 6 variables, Cb, A</?, 
Mp+£- , M e g, y/s miTl and Ef, all have a certain degree of sensitivity to spins, and are generally 
promising for spin studies. However, at this point one should not read too much into Fig. |2], 
since the observed variations in the kinematic shapes can be attributed to several different 
factors, not all of which are related to spins: 

• First, we see that in Fig. ^ there is a noticeable difference between lines of the same 
color, i.e. when the same spin scenario is shown for two different mass spectra. This 
means that the shapes of the kinematic distributions are affected by the mass spectrum 
(which in turn determines the available phase space) . This is not necessarily a problem 
for spin studies per se, since the mass spectrum will be known in advance. However, 
this does present a problem in the sense that any conclusions that one might draw as to 
which variables are most sensitive to spins, will necessarily be subject to the choice of 
mass spectrum. In other words, when it comes to spin determinations, it may very well 
be that for one particular mass spectrum a certain kinematic variable performs best, 
but for a different spectrum another variable is the winner. This is why in what follows 
we shall take great care in illustrating our results throughout the full mass parameter 
space, and not just a couple of study points as in Fig. || (see Sec. ^T] below). 



One should also keep in mind that the \fs distribution also varies from one model point 
to another, and this will impact the distributions of kinematic variables which carry a 
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Figure 2: Unit- normalized distributions of the nine variables introduced in Section |3| at a 7 TeV 
LHC. Each panel shows 4 cases, depending on the mass spectrum from Table |[ SUSY-like (dotted 
lines) or UED-like (solid lines); and the spin configuration: VSF for the MSSM (blue lines) and VFV 
for MUED (red lines). 



dependence on the underlying y/s in the event. Fig. || shows the y/s distributions for 
the four cases discussed in Fig. |2|. As expected, the SUSY spin scenario (VSF) leads to 
a significantly harder y/s distribution (blue lines in Fig. [|), due to the p-wave threshold 
suppression for scalar B1B2 production. Now returning to Fig. |2|, it should be rather 
intuitive that Cb, A</j, Mg+g-, M e g, \/s min an< ^ &t are au sensitive to the value of y/s, 
so it is not immediately clear how much of the differences between the red and blue 
lines in their distributions were simply due to the different y/s distributions in Fig. |3| 
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Figure 3: Unit-normalized y/s distributions for the four cases presented in Fig. |[ 

as opposed to spins. Conversely, the variables Mqt and Mqtx are less sensitive to y/s 
because of the back-to-back boost invariance. Therefore they will not be subject to the 
\fs effects in Fig. |3], and will have much smaller variations, as indeed observed in Fig. |2[ 

Of course, the shapes of the kinematic distributions do carry information about spins. 



As already explained in Sec. 2.2, this information can be extracted by comparing the 



predictions from the eight different spin scenarios in Table [l] to the observed data. 

• As emphasized in [20], spin correlations appear not only because of spins, but also due 
to the chirality of the fermion couplings. In reality, one is measuring the combined 
effects from spins and chirality, and the blueprint for disentangling these two effects 
can be found in [20,23,24,29]. We shall investigate the chirality effect in Section ||. 

4.1 Kinematic dependence on y/s and the mass spectrum 

We shall now discuss in some more detail to what extent our conclusions from the simple 
exercise in Fig. [2] persist as we vary the mass spectrum and the energy y/s in the event. 
In Fig. |I| we compare the VSF spin scenario of the MSSM to the VFV spin scenario in 
MUED, as a function of y/s (plotted on the x-axis) and the mass splitting between B and C, 
parameterized through the mass ratio Mc/Mb (plotted on the y-axis), for a fixed Mb = 500 
GeV. Each panel shows a temperature plot of the x 2 measure for one of the nine kinematic 
variables from Fig. [2|. Larger values of x 2 (warmer colors in Fig. ||) indicate larger differences 
in the predicted shapes between the SUSY and MUED scenarios, and smaller values of x 2 
indicate rather similar shapes, where the two scenarios might look the same. 

As already mentioned, the observed differences in shapes in general could be due to spin 
correlations, but could also be explained by different kinematics, as suggested by Fig. ||. In 
order to fairly assess the performance of the different kinematic variables with respect to pure 
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Figure 4: ^-comparisons of the two spin scenarios (VSF in the MSSM versus VFV in MUED), for 
each of the nine observables from Fig. |[ We fix the mass of B as Mb — 500 GeV, and then vary the 
CM energy y/s and the mass Mq of C. 



spin effects, one should compare the \ 2 values on the different panels in Fig. |] for a common 
value of y/s, thus taking kinematics out of the equation. Then, the variable with the larger 
X 2 for a given value of y/s would be more sensitive to spin correlations, and vice versa. 

There are several lessons to be learned from Fig. First, not surprisingly, different 
variables exhibit different degrees of sensitivity to spin correlations. In this sense, Mot, 
Mctx, Mt2 and Ei appear to be the worst performers. For the first three variables, this 
might have been expected, based on the result from Fig. |2[ The surprising member of this 
group is the lepton energy Ef. its distributions now appear identical in the two spin scenarios, 
in spite of the much more promising differences seen previously in Fig. |2[ This suggests that 
the differences in the E^ distributions in Fig. |2| can be attributed solely to the different y/s 
kinematics seen in Fig. |3| and not to true spin correlations. 

Fig. |] also reveals that the contrast between the two spin scenarios depends quite signif- 
icantly on the energy y/s, especially in the case of Cb, Aip, M e g and y/s min . In general, as 
the value of y/s increases, differences become more pronounced. For the largest values of y/s 
seen on the plots (several times above threshold) the most discriminating variable appears 
to be Cb, however such conclusion would be rather premature. As Fig. |3| shows, at hadron 
colliders particles are produced near threshold, with typical values of -^/i only 10-20% above 
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Figure 5: The same as Fig. [|, but instead of a linear collider at a fixed CM energy y/s, results are 
shown for the 8 TeV LHC. Here we vary both Mb and Mc and plot in the (Mb, -ff^-) plane. 



the 2Mb threshold. Unfortunately, Fig. |] shows that for those low values of y/s, the dis- 
criminating power of the Cb variable is significantly reduced, and it behaves similarly to the 
other alternatives like M e g and A<p. 

In order to get a sense of the full picture at the LHC (currently running at 8 TeV 
CM energy), one needs to re- weight the results from Fig. |] by the proper -y/i probability 
distributions analogous to Fig. ||| The result is shown in Fig. || where this time we vary both 
Mb and Mc and plot in the (Mb, j^) plane. Interestingly, Fig. |E| shows that in practice the 
variables whose distributions are most likely to show a noticeable difference between SUSY 
and MUED are M e g, \/s min and M^+£-. All three of these variables are sensitive to the 
yfs in the event, and are therefore capable of discriminating the different shapes of the y/s 
distributions in Fig. ||[ 



5. Expected shapes in the absence of spin correlations 



In this section, we shall be interested in analytical predictions of the shapes of some of the 
observables from the previous two sections. Consider a generic kinematic variable V. In 
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general, it is a function of the underlying event kinematics described in Sec. 2.1: 

v = v (VI, &*, ^ e[, ()>[■, e' 2 ,<f>' 2 ), (5.1) 

= V^*,®*,^,^,^). (5.2) 

Once this function is known, the pure "phase space" distribution of the variable V for a given 
fixed yfs can be obtained by integrating over dQ[ for the B\ — > t~C\ decay in the CMB1 
frame, over dVfl 2 for the B 2 —> i + C 2 decay in the CMB2 frame, and over dQ for the initial 
pp —> B\B 2 scattering in the CMBB frame: 

(V, yfs) oc J dtlj dfl[ J dfl 2 8 (V - V{^~s, &*, <p; 9[,<P[;6 2 , <j>' 2 )) . (5.3) 

Finally, at a hadron collider one also needs to convolute the fixed yfs distribution ( p.3[ ) with 
the parton luminosity L(yfs) (recall Fig. ||), so that the experimentally observed distribution 
will be given by 

diV f . _ , ^ f dN^ 



oc J dV~sL(V~s)(^)(V,V~s) (5.4) 

d^L(^) J dn J dn[J (m' 2 5 (v -v(y/s,e*,<p-, 9[, </>[-, ^fa)) . (5.5) 



dV 



Unfortunately, even if we concentrate on the fixed y/s case of ( [5. 3D , the 6 integrations in it 
are typically quite involved, and exact analytical expressions for the resulting distribution 
dN/dV{V,yfs) are known only in very few special cases (reviewed below). To make matters 



worse, the fixed yfs distribution of (5.3) still needs to be integrated numerically over the 
parton luminosities as in (|5.4|), Given all those difficulties, analytical studies of the kinematic 
shapes ( |5.5| ) (or even the simpler case of (|5.3|) ) appear to be quite challenging, especially for 
arbitrary values of yfs. 

In this and the following section, therefore, we shall limit ourselves to a much more 
manageable task: instead of dealing with a general yfs, we shall derive analytical formulas 
in some interesting and relevant yfs limits. There are two special values of yfs, namely 
the endpoints of its definition interval: the threshold value yfs t h = 2Mb and the infinite 
energy limit yfs —> oo. Taking either one of those limits leads to simplifications in the 
defining functions ( |5.1| , ^2] ), and the integrations become easier, although success is not always 
guaranteed. We shall now discuss selected kinematic observables for general yfs (where 
possible), and in the two limits of yfs —> xfs^ and yfs —> oo. We have already learned from 
Fig. H that out of those two limiting scenarios, the threshold case is much more relevant at 
the LHC, since the mother particles B are much more likely to be produced near threshold 
as opposed to very large boosts. 



5.1 The relative pseudorapidity variable cos(9 



We begin with the Cb variable from Section 3.1 



C B = tanh ( ^ — ^ ) . (5.6) 
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The lepton pseudorapidities in the LAB frame 771 and 772 (which were defined in ( |3.3D ) can 
be conveniently expressed in terms of the respective pseudorapidities ( |2.8|) in the CMB1 and 
CMB2 frames as follows 

(— 1)*77*) + cos©*sinh (77- - (— l)*r/*) + cos sin©* _ _^ 

N sinh (77' — (— 1)'77*) — cos <bt- sin 0* 



1 



Vi = Pz + 2 ln 



cosh (77^ - v ^ ., y , ^ „ ^ v ., t v ^ ., L , ^ 
cosh (77^ — (— 1)*7/*) — cos 0* sinh (77^ — (—!)*?/*) — cos 1 



Substituting (|5,6|) and (5/7) into (|5.3j ), one can in principle obtain the Cb distribution for a 
given fixed y/s. Unfortunately, we have not been able to obtain closed form expressions for 
general y/s, and shall consider the two limiting cases instead. 

In the infinite energy limit y/s — > 00, eq. (|2.6| ) implies that 77* — > 00 and then (|5~7i 
reduces to 



m 



-> Pz + In 



-> p z + In 







cot 1 












tan 









Then (5.6) simplifies to 



C B 



-> cos ©* 



(5. 



(5.9) 



Therefore, in the infinite energy limit, the Cb distribution reduces to the cos 0* distribution: 

dN dN 



dCB ??*^oo dcos0* 



(5.10) 



This fact served as the original motivation for introducing the Cb variable in the first place [8] . 
It is well known (see Section |6.2| below) that the cos 0* distribution is directly probing the 
spins of the particles A and B, and so Cb will inherit some spin sensitivity through ( 5.10| ). 

The correlation (|5.9| ) is pictorially illustrated in Fig. ||, which shows scatter plots of Cb 
versus cos 0* for several different boost factors 77* for the two study points from Table g. We 
see that ( jsT9| ) holds only as long as y/s is sufficiently large. As y/s gets closer to threshold, 
the correlation is lost and so is the sensitivity of the Cb variable to the spin effects encoded 
in the cos©* distribution (this effect was already evident in Fig. ||). 

Given that most events at a hadron collider are produced at threshold, we now turn to 
the opposite limit of y/s ~ \/s th or equivalently, 77* — > 0. In that case, the angle 0* becomes 
arbitrary and can be set to zero, so that ([5l]) and ( |5.7| ) give 



Cb 



>o,e*->o 



-> tanh 



V'2 



(5.11) 



With this simple result, one can now perform the integrations in ( |5.3| ) and obtain the unit- 
normalized Cb distribution at threshold as 



dN 



1 



C 2 



d Cb v* 



4 CI 



-2C B + (1+ Cl) ln 



1 + C B 
1- C B 



(5.12) 
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Figure 6: Correlation between Cb and cosO*, for different boost factors r/*. The upper (lower) row 
of plots is for the MSSM (MUED) study point from Table § 

This formula is one of the main new results in this paper. It provides the benchmark shape 



of the Cb distribution in the vicinity of the relevant energy regime. Eq. ( |5.12D will be central 
to our understanding of the effects of spin correlations later on in Section ||. Of course, 



( 5.12 ) is nothing but a crude approximation to the true Cb distribution which is obtained 



by integrating over all energies yfs as in (|5.5|) . Nevertheless, for practical purposes ( 5.12 



appears to be a more relevant limit than fl5,10|) , since the typical values of yfs are closer to 
threshold than to oo — see Fig. ||. 

5.2 Azimuthal angular difference Aip 

The azimuthal angular difference A93 defined in Sec. |3.2| is given by 



where 



cos A0' + A) 

cosA(£ = — — , (5.13) 

L)\L)2 



Dq = sin 2 O* (cos cos (j)' 2 — sinh r][ sinh rf 2 
1 

~ 2 



- sin 20* (cos 4>2 sinh r( x — cos 0j sinh rj' 2 ) , (5-14) 



Di = {1 + sin 2 9* (sinh 2 f}[ - cos 2 $[) - sin2G* cos^ sinh f)[} 1/2 , (5.15) 
D 2 = {1 + sin 2 9* (sinh 2 fj' 2 - cos 2 cf) 2 ) + sin2G* cos^ 2 sinh^} 1/2 , (5.16) 
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and 



77 +77!, 

Now in the infinite energy limit we obtain 

Aip 



and the universal result 



diV 



dA<p r)* 

In the other y/s extreme, at threshold, we obtain 



fj' 2 =r,*- V f 2 . (5.17) 
-^vr (5.18) 

>5(A<p-ir). (5.19) 
n 

> Acp' (5.20) 

const. (5.21) 

dA<£> »?*->0 7T 

We see that as y/s varies between its two extremes, the Aip distribution changes from being 
completely flat to a perfect delta-function peak at tt. Therefore, the flatness of the A(p 
distribution is a measure of the typical y/s in the event, which in turn is affected by spins as 
shown in Fig. [| Therefore one might expect that the A(p distribution would also be sensitive 
to spin correlations through y/s effects, and this is indeed what we saw previously in Figs. |2[ 
| and |. 

5.3 Invariant mass Mg+ e - 

On our notation, the formula for the rescaled invariant mass variable (3.7) is 



and a flat distribution 



dN 



M l+l - 



I cosh (r)[ + r]' 2 ) — cos (4>[ 
2 cosh ?][ cosh 772 



(5.22) 



The integrations in (|5.3|) can be performed for any fixed y/s > 2Mb and the result for the 
unit-normalized M e +£~ distribution is [36] 



diV 



dM en - 



sinh 1r\ 



-M e+e -, 



sinh 2»7 



-M, 



i+t 



log 



M T+t- 



M 



where 



M cusp 

lvl e+£- 



for M i+i - <M™/_, 



for M^ p < M t 



M £+e- 



< A/f max 



(5.23) 



(5.24) 



For events with a given fixed value of y/s, the distribution ( 5.23| ) exhibits a "cusp" , i.e. a non- 
differentiable point, at Mg+£~ = M™ S F_ [36]. However, since the events at hadron colliders 
occur at a variety of values of y/s (see Fig. |3|), the cusp gets smeared, since its location 
depends on the boost factor 77* and thus on y/s. 
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Figure 7: Top: \ 2 tests of the predicted Mctx shape (5.25) against simulated data in MUED (left) 
and SUSY (right) for a fixed M B = 500 GcV and varying </s and M C /M B - Bottom: the M C t x 
distributions for the study point So with the largest % 2 . 



5.4 Doubly projected contransverse mass Mqtx 

In the absence of spin correlations, the analytical formula for the Mqtx distribution (re- 
stricted to events with non- vanishing Mqtx) is given by [41] 

dN 4M CT ,mfeV (,25) 



dM CTx \M ( ™ x x \ 

Interestingly, this formula does not carry any explicit dependence on y/s, which can be 
traced back to the invariance of the Mqtx variable under back-to-back transverse boosts 
of the parents B\ and B2 [41,55]. Even more peculiar is the observation made in Figs. |2[ 
|I] and [|: that the shape of the Mqtx distribution appears to be pretty much independent 
of the spins and the mass spectrum in the model. One may therefore wonder whether the 
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distribution ([5.25 ) is indeed a universal function which is completely independent of the spin 
and the kinematics. 

In order to address this question, in Fig. [?] we test how well the formula ( 5.25| ) fits 
the Mqtx distributions in the MUED and SUSY spin scenarios, over a wide range of mass 
spectra. The upper two panels perform a \ 2 test of the predicted shape ( p.25[ ) against the 
Mqtx distribution obtained in our simulations for MUED (left panel) and SUSY (right panel) . 
The test is repeatedly performed for various fixed values of y/s (plotted on the x-axis) and 
mass hierarchies Mq/Mb (plotted on the y-axis), for a fixed Mb = 500 GeV. We see that in 
the MSSM, the Mctx shape is very well predicted by Q5.25 ) in all cases. This could have been 
expected because of the scalar nature of the intermediate particles B, which causes the two 



leptons to be uncorrelated. The formula (5.25) works very well in MUED as well, and only 
for Mc « Mb one starts to see deviations. In order to illustrate the size of the deviations, 
in the lower two panels of Fig. ^ we plot the Mqtx distributions for the point So with the 
largest x 2 in the MUED case: Mc = 0.1Mb and y/s = 14.4Mb- We see that even in the 
worst case scenario of So, the two distributions are in reasonable agreement. We conclude 
that (5.25) is not a universal function and the shape of the Mqtx variable in principle does 
get affected by spin correlations, although the effects are very minor. 

5.5 Lepton energy Eg 

The distribution of the lepton energy Eg has the famous box-like shape for any given fixed 



diV 
dE e 



l 

0, 



if E { r n) <E f <M max \ 



otherwise, 



(5.26) 



where Eg is the rescaled energy variable from (3.17). The two endpoints of the distribution 

(5.27) 



(5.26) are given by 



E 



{rain) 



E 



(max) 



In the threshold limit of y/s — > 2Mb or rf — > 0, ( |5.26 ) becomes simply 

diV 



dE f v*^o 



> S(E t - 1). 



(5.28) 



The observation that the two lepton energies are equal and constant at threshold was the 
main inspiration for the razor class of variables [61]. At hadron colliders, where y/s is varying 
from one event to another, the distribution ( |5.26D needs to be convoluted as in ( |5.5D and the 



two endpoints (|5.27|) become smeared, as can be seen in Fig. g. 

The distributions of the remaining four variables from Sec. |3] (Mcr, M e g, y/s min and 
Mt2) do not have compact analytical expressions and we shall not discuss them here. 



6. Spin correlations in the Cb distribution 

In this section, we shall focus on the Cb variable and investigate how spin correlations affect 



the shape of its distribution. In Section 5.1 we saw that even in the pure phase space 
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approximation (with no spin effects), we managed to obtain analytical formulas in closed 
form only for the two limiting cases of \fs at threshold and at infinity. Correspondingly, we 
shall now study the effect of spin correlations on Cb in those two limits as well. 

6.1 Spin effects at threshold 

We begin with the threshold limit y/s ~ (= 2Mb), where the Cb distribution is given by 
the previously derived formula ( |5.12| ). Then we ask, how will the presence of spin correlations 
modify the phase space result ( |5,12| ). For this purpose, in Appendix ^ we rederive the unit- 
normalized Cb distributions at threshold for each of the 8 spin scenarios in Table [l], and for 
arbitrary values of the relative chirality parameters a, /3 and 7 defined in ( [2.12f^T4|) . The 
result is 
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where 



Ml 
Ml 



and the basis functions Jj are defined as follows, 
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Figure 8: The basis functions J(Cb) denned in (6.10-6.14). The solid red line is the pure phase space 
distribution ( 6.1 0| ) or equivalently ( |5. 12 ). 
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(6.13) 
(6.14) 



The basis functions ( |6.10| - 6.14 ) are pictorially illustrated in Fig. ||. The function Jps is simply 
the phase space limit (|5.12j ) already derived in Sec. 5.1 , while the remaining 4 basis functions 
parameterize the different possible distortions of the pure phase space shape, for the various 
spin scenarios. 

Formulas (6.1-3.8) represent one of our main results. Using fl6.lH6.8p , one can understand 
how the shape of the C B distribution near threshold changes as a function of the mass 
spectrum (through the dependence on the y parameter) and as a function of the relative 
chirality parameter 7 defined in ( 2. 12] ). (The parameters a and /3 decouple and do not enter 
the threshold distribution formulas.) For example, ( |6.1| ) and fl6.5| ) show that in the cases 
where the intermediate particles B are scalars, any spin correlations between the two leptons 
are wiped out and one obtains the pure phase space shape ( |6.10| ), which peaks at C B = 
and vanishes at C B = ±1 (see the solid red line in Fig. ||). Furthermore, the predictions ( |6.1| ) 
and ([0]) are not sensitive to the mass spectrum at all, since the mass ratio y does not enter 
those formulas. 
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Figure 9: Spin correlation effects on the Cb distribution at threshold, for various spin scenarios. 
In each case, we start with the pure phase space distribution at 7 = and/or y = 1, and show the 
amount of shape distortion by varying 7 and y over their full range. 



Now let us proceed to the SFS and VFS scenarios of ( p.2|) and (p.q), correspondingly. 
Here we find that the spin correlations induce a second term involving the Jf function, which 
has a maximum at Cb = and two local minima at Cb = ±0.845 (see the black dashed line 



in Fig. g). The coefficient of the Jf term is positive in (6.2) and negative in ( |6.6D . This means 
that in the SFS spin scenario, the addition of spin correlations will make the central peak of 
the pure phase space distribution Jpg even more pronounced. This is illustrated in the upper 
left panel of Fig. [9|, where we start with the pure phase space distribution at 7 = (the 
black solid line) and gradually increase the amount of spin correlations until we get purely 
chiral couplings at 7 = 1. As expected, the presence of spin correlations in the SFS scenario 
makes the central peak in the Cb distribution steeper. In the VFS scenario, on the other 
hand, spin correlations have just the opposite effect, since the Jf term in ( |6.6D comes with a 
minus sign. Now, the central peak will be suppressed, and events will shift towards the two 
endpoints instead. The resulting shape distortion is illustrated in the upper right panel of 
Fig. [|, where again we show the full range from 7 = (pure phase space) to 7 = 1 (purely 
chiral couplings). We see that the presence of chiral couplings causes the VFS distribution 
to be much flatter than in the pure phase space limit. 

The SFS and VFS results shown in the upper two panels of Fig. |9| are also independent 
of the mass spectrum, as the mass parameter y does not enter eqs. (|6.2|) and (|6.6[). In 



contrast, the remaining four spin cases to be discussed now will exhibit a dependence on 
the mass spectrum through y. It is probably easiest to start with the SFV and VFV spin 
configurations, whose Cb distributions are given by ( |6.3[ ) and (|S7^), respectively. We find very 
similar behavior to the two cases of SFS and VFS just discussed — the only difference is that 
the coefficient of the Jf function is now additionally suppressed by a factor of (1 — 2y) 2 /(l + 
2y) 2 , which incorporates the mass spectrum dependence. Because of this suppression, one 
would generally expect smaller variations than what we saw in the upper two panels of 
Fig. ||. There even exists a very special case with y = 1/2 (i.e. Mb = y/2Mc) when the 
Cb distribution reduces to the pure phase space prediction ( |6.10|) . For that particular mass 
spectrum, the SFV and VFV distributions at threshold would be identical to the SSF and 
VSF distributions, and spin discrimination would have to be done in a regime away from 
threshold. 

Let us now move on to the SVF scenario of ( |6.4D . This is perhaps the most complicated 
case, since in addition to the Jps baseline shape it has two extra terms — one with the 
already familiar function Jf, and a second term involving the function Jy (shown in Fig. |8| 
with the blue dashed line). The function Jy peaks at Cb = and has broad shoulders near 
the endpoints of the Cb interval. The coefficient of Jy in fl6.4| ) is always positive, which 
means that the Jy term will always tend to enhance the central peak. At the same time, 
the coefficient of the Jf term in (6.4) is always negative, which would tend to suppress the 
central peak, as already seen in the upper right panel of Fig. ||. Thus the interplay of the two 
terms will lead to partial cancellations, which are illustrated in the lower left panel of Fig. |9]. 
The dashed blue line represents the pure phase space limit in eq. (|6.4D , which is recovered 
when 7 = and y = 1. If we keep 7 = 0, the Jf term is turned off, and by varying y, we 
can see the impact of the Jy term in ( |6.4[) , which is maximal at y = (the solid blue line). 
Conversely, if we keep y = 1, the Jy term is turned off, and the variations in 7 reveal the 
effect of the Jf term, which is maximal at 7 = 1 (the red dashed line). The case when both 
terms are maximal (7 = 1 and y = 0) is given by the solid black line. 

Finally, we discuss the VVF case of eq. (|6.8| ), pictured in the lower right panel of Fig. |9[ 
Here the Cb distribution at threshold depends only on the mass spectrum but not on the 
chirality parameter 7. The spin correlation effects are encoded in the two functions Jyi and 
Jy2, represented in Fig. || by the green dotted and magenta dot-dashed lines, respectively. 
Jyi has local maxima at Cb = and Cb = ±0.881 and local minima at Cb = ±0.456, while 
Jy2 has local maxima at Cb = ±0.7 and a local minimum at Cb = 0. The lower right panel in 
Fig. |9] shows the resulting variation in the Cb shape as we vary y from y = 1 (the pure phase 
space limit) to y = (when the Jyi term is maximal). We see only very minor variations, 
due to the fact that the coefficients of both the Jyi and Jy2 terms in ( |6.8[ ) are numerically 
very small. 

This concludes our discussion of eqs. (0-3^), which manifestly describe the spin ef- 
fects at threshold. As Fig. || showed, although most events are produced near threshold, 



virtually none are produced exactly at threshold, thus eqs. (6.1-|6.8j) are perhaps of limited 



practical interest. Their true value is in developing some intuition about the shape of the Cb 
distribution, which could prove useful for the interpretation of the finite ^fs results below. 
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6.2 Spin effects in the large energy limit 



We now turn our attention to the other y/s extreme, namely y/s — > oo. In that case, as dis- 
cussed in Sec. 51 and shown in ( 5.1 0| ) , the Cb distribution reduces to the cos0* distribution 
for 2—7-2 processes, whose dependence on the spins of the particles is well known, so we have 



diV 



dC B 



SXY 



diV VFY 

d Cb oo 
dN VYF 



dC B 



oc 1, for {X,Y} C {S, F, V} 
oc 1 + C B , for Y G {S, V} , 
ocl - C B , for Y G {S, V} . 



(6.15) 
(6.16) 
(6.17) 



These equations predict three generic shapes for the Cb distribution at large energies: 



Flat as in ( 6.15 ). This is the case whenever A is a scalar particle. 



1 + C| as in ( ISleD - This case occurs when A has spin 1 and B is a fermion, e.g. the 
UED example from Section |3] falls into this category. The distribution ( |6.16| ) peaks at 
Cb = =tl and has a minimum at Cb = 0. 



1 — Cg as in fl6.17p . This is the case when A has spin 1 and B is a boson, as in the 
SUSY example considered earlier. The Cb distribution then has a peak at Cb = and 
vanishes at the endpoints Cb = ±1. 



As shown in Sec. 5.1, the nice correlation ( |5.10| ) between Cb and cosB* emerges only at 
sufficiently large energies, which might not be realistically achieved at a hadron collider. 
Therefore, the practical value of the limits ( |6.15f]6~T7| ) is also debatable, just like the formulas 
( 6.1 -|B78 ) from Sec. |B~l] . However, with the use of ( |B~l] - |6.8D and ( |6. 15 -6. 17), one can begin to 
understand the observed Cb shapes in Fig. ||, where we noticed that the MUED case gives 
a flatter Cb distribution than the MSSM. In the case of the MSSM, the Cb distribution at 
threshold (|6.5[) peaks at Cb = and vanishes at Cb = ±1, while the large energy limit (|6l7D 
also peaks at Cb = and vanishes at Cb = il. The Cb distribution at finite y/s would then 
interpolate between those two very similar limiting cases and inherit their common properties 
— which explains why the actual MSSM distribution in Fig. [2| has a relatively sharp peak at 
Cb = and vanishes at Cb = ±1. 

In the case of MUED, we have a different story — here at intermediate y/s we have to 
interpolate between two very different shapes: the threshold distribution fl6.5|) which peaks 
at Cb = and vanishes at Cb = ±1, and the asymptotic distribution ( |6.16 ) which peaks 
at Cb = ±1 and has a minimum at Cb = 0. At the end of the day, the resulting Cb shape 
will depend on the exact form of the underlying y/s distribution, but Fig. [3| suggests that the 
threshold behavior will dominate, since a typical event is produced closer to threshold. This 
is indeed what we observe in Fig. [2], where the Cb distribution in the MUED case also peaks 
at Cb = 0, in spite of the asymptotic expectation ( |6.16 ). 
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6.3 Spin effects near threshold 

When particle A is off-shell (i.e. Ma < 2Mb), the onset of the pp — > B\B\ production cross- 
section as a function of yfs is sensitive to the spins of A and B. The partonic cross section 
at a given ^fs can be expressed as 

H^) = yJ dn 2l M (^)! 2= (i) 2 ^f / d^|M(v^)| 2 , (6.18) 

where /3cm is the B{ boost factor in the CMBB frame 

/3 CM ^tanh??*. (6.19) 

The cross-section near threshold is suppressed by a certain power of /3cm, which depends on 
the spin scenario: 



\ M (^)\%s) °cl + 0(/3 CM ), ( 6 - 2 °) 
|M(v^)|( SF) «(l-cosJ) + (l + co S( 5)/3 CM + 0(/3 CM ), (6.21) 

l M (^)l?sv) °c1 + 0(/3 C m), ( 6 - 22 ) 
|M(Vi)|( VS) «/3 C M(l-cos 2 e*)+0(/3 CM ), (6.23) 

l M (v^)|( VF) ocl + 0(/3 CM ), (6.24) 
I M (^)I(W) K ^cm (l - ^ cos20 *) + °(^cm) • (6-25) 
Interestingly, the case can be subdivided into two categories: a real scalar 

|M(VS) |( S++F) oc /3 CM + 0(/3 CM ) , (6.26) 

or a pseudo-scalar 

|M(^)|J s _ +f) oc1 + 0(^cm)- (6-27) 
Thus, the /3cm suppression factor near threshold is as follows: 

(6 - 28) 

[/3 CM , for(S^F) = (S ++ F),(VS),(VV). 
7. Comparison of the different spin scenarios 

We are now ready to contrast the different spin scenarios and discuss the prospects for spin 
discrimination. In principle, one could pose two questions: 

1. How well can two different spin scenarios be discriminated experimentally, given the 
uncertainty in the coupling chiralities and instrumental effects like SM backgrounds 
and the finite detector resolution. Clearly, the answer to this question will depend on 
many quantitative factors - the chosen study point, the size of the signal, etc., and it 
is difficult if not impossible to give a"one size fits all" answer. 
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VFS 


2 

Jps — -4- Jf 




1 


VFV 


VFV 


2 

Jps — -4- ei Jf 


i+el 


1 


VFS 



Table 4: The asymptotic behavior of the Cb distribution at threshold (second column) and at 
yfs — > 00 (third column) for different spin configurations in the antler topology. The third column 
lists the power of any additional /3cm threshold suppression coming from the matrix element. The 
last column shows the alternative spin scenario with similar properties. 



2. Which spin scenarios (and under what circumstances) are in danger of being confused 
with each other? This question is easier to tackle theoretically, because if we can identify 
the cases where two spin scenarios look identical at the parton level, the conclusions 
will remain unchanged when we add all the usual experimental complications. This is 
why in this section we shall use the intuition developed in previous sections to pinpoint 
the difficult cases for spin discrimination. 

Let us focus on the Cb distribution. Table § summarizes its salient features identified 
in Sec. [| the behavior at threshold (second column) or asymptotically at yfs — > 00 (third 
column), and the power of any additional /?cm threshold suppression arising from the matrix 
element. For convenience, formulas (^T-^8) are rewritten in terms of the shorthand notation 



ei -virW £2 "4(TT^y f3 -(iw (7 - 1} 

_ (1 - y) 2 y(l-y) 

64 " WTV) 65 " W+y)' (7 - 2) 

In Table ||, the 8 spin configurations from Table [l] are arranged in pairs. The two spin 
scenarios within each pair can exhibit identical behavior both at threshold (for a suitably 
chosen mass spectrum) and at yfs — > 00, and have the same /3cm threshold suppression. 
Those pairs are therefore the problematic cases where spin discrimination based on Cb might 
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Figure 10: The same as Fig. ^, but this time comparing kinematic distributions for the SSF and SVF 
spin scenarios at an 8 TeV LHC. The couplings and masses in each case are the same, and were chosen 
to match the asymptotic Cb distribution at threshold (y ~ 1 and 7 = 0): Ma = 1200 GeV, Mb — 500 
GeV and Mc = 495 GeV. For this case, (3 is not defined, while the choice of a is inconsequential. 



be difficult. Of course, this is simply a conjecture based on the asymptotic behavior, so it is 
worth checking if it holds for the Cb distribution at any y/s. Then if it turns out that the Cb 
distribution is indeed very similar for a given pair of spin configurations, we will then check 
whether any of the remaining variables discussed in Sec. || can offer an alternative tool for 
discrimination. 

Figs. [lO|-|r^ show the distributions of the nine variables from Sec. ||[ for each of the twin 
spin scenarios in Table ||[ All simulations were done with MadGraph v4 [62] , in the narrow 
width approximation, for an 8 TeV LHC. In all figures, we consider the case of a heavy 
resonance A with Ma = 1.2 TeV, and intermediate particles with mass Mb = 500 GeV. 
We then adjust the mass of Mc appropriately in order to obtain an exact match in the Cb 
distributions at threshold. For example, Fig. 10 compares the case of SSF to SVF. The former 
has a Cb distribution which is given by the pure phase space formula at threshold, while the 
Cb distribution for the latter involves terms proportional to €2 and e37 2 /4 (see Table 0). 
Then the worst case scenario would have a degenerate spectrum with y — > 1 (Mg ~ Mc), 
which would set €2 = 0, and purely vector-like couplings with 7 = 0. This is precisely 
the case shown in Fig. 10, where for definiteness we choose Mc = 495 GeV. As expected, 
the two Cb distributions are almost identical and can hardly be used to measure the spins. 
The remaining plots in Fig. 10 then check whether any of the other eight distributions show 
observable differences. As it turns out, they are all pretty well matched as well, and the spin 
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Figure 11: The same as Fig. |lCj, but for the SFS and SFV spin scenarios with a scalar particle A. 
Here the threshold behavior for Cb is obtained with y ~ 0, so unlike Fig. [l(], here we choose Mc = 5 
GeV. 



discrimination between SSF and SVF is indeed very problematic under these circumstances. 
In the next two figures we show a similar comparison between the SFS and SFV spin 



scenarios: in Fig. 11 for the case when particle A is a scalar (S ++ ) and in Fig. [T^ for the case 



when particle A is a pseudoscalar (S + ) Here, to match the threshold Cb distributions, it is 



sufficient to set y = (for definiteness, we choose Mc = 5 GeV). Figs. 11 and 12 reveal a 
similarly grim situation — the distributions of all 9 variables match almost perfectly in the 
two cases, and any spin discrimination thus appears to be virtually impossible. 

The next pair of similar spin configurations in Table |I] is VSF and VVF, which are 
compared in Fig. The threshold Cb distribution for VSF is given by the pure phase space 



result ( 6.10| ) while for VVF, it reduces to the same formula if y = 1. In order to mimic the 



case of y = 1, we again choose Mc = 495 GeV. Once again, we find that the distributions of 
all 9 variables are pretty similar. 

The last problematic pair of spin configurations in Table f| is VFS and VFV. Their Cb 



distributions at threshold become identical if 7 = or if y ~ 0. Fig. |FJ shows a case with 
7 = and y ~ 0. Once again we observe that all 9 kinematic variables behave very similarly 
and spin discrimination is very difficult. 

In conclusion of this section, we can summarize its main two lessons as follows: 

• For the spin scenarios where particle B is a fermion, the two situations in which it will 
be quite difficult to determine the spin of particle C are: 



- 33 - 





Figure 12: The same as Fig. [Tl], but for a pseudoscalar particle A. 

1. Vectorlike couplings between B and C with 7 = 0. The presence of chiral couplings 
is a necessary condition to observe spin correlations [12,18,20,27]. 

2. Models with y ~ 0, where particle C is highly boosted. When C is a vector boson, 
it is longitudinally polarized and effectively behaves as a scalar. 

• For the spin scenarios where particle B is a boson, it may be difficult to determine 
whether that boson is a spin or spin 1 particle. Models with 7 = and/or y ~ 1 are 
again the most problematic. 

8. Discussion 

In this paper we investigated the antler event topology of Fig. [l] with regards to spin effects 
in the kinematic distributions of the two visible particles. Given the simplicity of Fig. [I], one 
might have hoped to be able to describe analytically the observable kinematic distributions, 
including explicitly the effects of spin correlations. Unfortunately, this is not the case, and 
spin studies must be done numerically by template methods. Nevertheless, useful physics 
intuition can be developed in the two -y/s extremes, namely, at threshold -y/s ~ 2Mb and 
at -^/s — > 00. In Section || we analyzed the asymptotic analytical expressions for the Cb 
distribution, which helped us identify in Sec. the difficult cases for spin discrimination. In 
Table |3] we identified pairs of spin configurations which can easily mimic each other, and 
the exact circumstances when this is most likely to occur. Typically, the chirality 7 of the 
BCi coupling and the mass splitting between B and C play an important role — the spin 
correlations tend to be diminished when 7 = and depending on case, 
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Figure 13: The same as Fig. |Io|, but for the VSF and VVF spin scenarios with y ~ 1 (Mc = 495 
GeV). 

• For the spin scenarios where particle B is a fermion: Mc ~ Mb- 

• For the spin scenarios where particle B is a boson: Mc « Mb- 

While our discussion centered mostly on the Cb variable, which has been most often 
utilized for spin studies, we also considered 8 other kinematic variables, which were previously 
developed mainly in the context of mass determination. Using the popular examples of SUSY 
and MUED, we investigated their power for spin discrimination. We generally found that 
some variables are better suited for spin studies than others. Among the better performers 
were Cb, Ay?, Mg+g-, M c g and y^mm- ^ ^ s interesting to note the complementarity between 
the Cb variable and the others in this group. As shown in Fig. ||, the Cb distribution becomes 
most sensitive to spins at high energies, where the correlation ( |5.10| ) is manifest. In contrast, 
the distributions of the other variables are probing the different ^/s behavior near threshold 
(see, e.g. Fig. |3|). Therefore it appears that a multifaceted approach, utilizing a number of 
different and complementary variables, would be most beneficial for spin determinations. 
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Figure 14: The same as Fig. [To[ but for the VFS and VFV spin scenarios with a = ft = 7 = and 
y ~ (M c = 5 GcV). 

A. Derivation of the Cb angular distribution at threshold 

In this Appendix, we show the derivation of the Cb angular distributions in the threshold 
limit \fs — > 2Mb- According to ( |5.11 ), in that case the Cb variable can be expressed in terms 



of the pseudorapidities of the visible particles in the corresponding CMB1 and CMB2 frames 



as 



C B =tanh(2i_!fc), (A.l) 



with and rj' 2 defined in (|2.8| ). Since the initial state at the LHC is symmetric (both beams 
are proton beams), one should symmetrize the Cb distribution with respect to Cb O — Cb, 
since we do not know which proton beam the initial state quark came from. 

A.l Cb distribution with no spin correlation 

In the pure phase space limit, 
dN 



(W oc JJJ dO^dOid^lMp^CB-CBfe*,^,^,^,^]) 

oc j j Arii dr?2 sech 2 r\ x sech 2 rj' 2 5 f Cb — tanh - ^ 2 ^ 

oc I dr/2 {cosh ^ (0) - 7/2) + l} sech 2 rj'^ sech 2 rj' 2 , (A.2) 
J— 00 
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where rj'i' is the zero of the argument of the delta function 



Cb — tanh 



Vi 



I (o) 



% +ln 



1+ Cb 
1- C b 



(A.3) 



After integrating over 772, we get the following unit-normalized distribution 



Jps 



dN _ 1 - Cg 



dC B " 4Cg 



1 + C B 



1- C B 



(A.4) 



-2C B + (1+ Cl) In 

which is eq. ( 5.12| ) . This result is already symmetrized with respect to Cb — Cb • 

A. 2 Cb distribution when particle i? is a fermion 

When particle B is a fermion, the matrix element has the generic form 

\M\ 2 oc 1+ A (cos 0i -cos^ + ^icosA^' + ^sinA^sin^sin^+Ccos^cos^. (A.5) 

The B\ and B2 terms integrate out to zero and we get 
dN 



dC B 



oc 



dfTd^df^ \M\ 2 5{C B -C B [e*,<^,7f,4^]) 



1 ( dTVps , diV^ r dAc\ 

iv tot v dc B dc B dc B y ' 



where 



dAps 
dC B 



oc^y dn[dn' 2 5(C B - C B [r&,r&]) 



(A.6) 
(A.7) 



is the pure phase space contribution already derived in Sec. A.2| , while the remaining two 
contributions are 



dN A 
dC B 

dNc 
dC B 



oc J J dOj d$7' 2 (cos 6'i — cos 6*2) 5 ( Cb — Cb [t/i, r) 2 ] ) > 
oc dO; dO' 2 cos 0; cos 6' 2 5 ( C B - C B [ryi , r/ 2 ] ) . 



(A.8) 
(A.9) 



The contribution from (A.8) vanishes after the symmetrization Cb — Cb and we are left 
to evaluate 



dN c 
dC B 



OC 



OC 



J J dCl[ dCl' 2 cos cos 9' 2 5 ( C B - C B f? 2 ] ) 

dr/i dr/2 sech 2 77^ sech 2 r/ 2 tanh rj'i tanh r? 2 5 ( Cb — tanh 



7][ - rf 2 



dr]' 2 |cosh (jj'i® — r/ 2 ^ + l| sech 2 rj[^ tanh 77^°^ sech 2 rj' 2 tanhr/ 2 



1 - C 2 , 

1 "^b J — 00 



sech 2 77 2 tanh r]' 2 tanh ^r/ 2 + In 


"l+C B " 
_1-C B 


) 


cosh 2 (^r] 2 + In 


_1-C B _ 


\ 


) 



(A.10) 
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Upon integrating over 773, this gives the contribution which is due to spin correlations, we 
use a subscript F for it as a reminder that it corresponds to the case of a fermion particle B: 



dN 1 fdN P s , „dA F \ C 1 



(A.ll) 



where ( [A.10| ) evaluates to 



dN F 1 - C\ 



dC B ACl 



6C B + 4C B + 6C B -(3+ C B + C B + 3C B )ln 



1+ C B 
1- Cb 



, (A.12) 



which is the result in (6.11). The prefactor C in (A.ll) depends on the particular spin 
configuration, as seen in eqs. (|Q,3.3, |6^J6.7| ). 

A. 3 Cb distribution in the SVF scenario 

In the SVF spin scenario, the matrix element has the form 

(i - y? 



\M\ l oc 1 + 



+ 



l + 2y 

27 2 



(cos A(f>' sin 0[ sin 6' 2 — cos 9[ cos O^Y 



l + 2y 



cos A(j)' sin 9[ sin 02 — cos 9[ cos 2 ) • 



(A.13) 



The Cb distribution is 
diV 



dC B 



oc 




1 /diVpg (l-y) 2 dA v 2 7 2 dA F \ 
A^ V dC B + l + 2y dC B ~ l + 2y dC B ) ' 



(A.14) 



where dC ^ s is the pure phase space contribution ( |A.4j ), while is the contribution ( A.12j ) 
from the helicity structure of the BtC interaction. The new term here is given by 

oc JJ dfiidfi^ (cosA^sin^sin0 2 -cos^cos0 2 ) 2 5 ( C B - C B [r/i, r/ 2 ] ) . (A.15) 

After integrating out we get 

oc J J d cos 0j d cos 2 Q sin 2 6[ sin 2 2 + cos 2 6' x cos 2 2 ^ 5 ( C B - C B [r?i , r/ 2 ] ) 



OC 



dr/i dr/ 2 sech 4 ry^ sech 4 ?? 2 ( — h sinh" 4 77^ sinh^ r]' 2 I <5 ( Cb — tanh 



1 



2 „/ „;„l.2 / 



drj' 2 I cosh (rj'i® — rj'^j + l| sech 4 sech 4 r/ 2 ^- + sinh 2 r/^ ' 1 sinh 2 r/^ 



1 - C 2 , 

^ "-B J —00 



sech 4 ?? 2 ^5 + sinh 


2 r/ 2 sinh 2 ^ 2 + In 


"l+C B " 
1-C B 
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cosh 4 


rfc + In 


_1-C B _ 







(A.16) 
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Finally, after integrating over r)' 2 , we get ( |6.4| ) 

dN _ 1 / diVps (l-j/) 2 dAT v 2 7 2 djV F \ 
~dC^~N^ t V dC B + l + 2y dC B ~ l + 2y dC B J 



Ml±^)( Jps + iLzC Jv 



7 



(2 + y) z V 4(1 + 2!/) 2(1 + 2y) 

where the new contribution is given by ( |6.12| ) 



Jf , 



(A.17) 



Jv 



dNy _ 1 - C\ 
dC B ~ 48 Cl 



2C B (45 + 22C B + 45C^) 



+ 3(l5-5Cl+6C B + 6C B -5C| + 15C B °)ln 



1+ C B 
1- C B 



(A.18) 



A. 4 C B distribution in the VVF scenario 

After integrating out from the matrix element, we get 

d^d^lMp oc 1 



U -y) 2 „„„2 fl / „„„2 fl / vO--y) i .2,,/ 



i + y 



COS COS P 2 



2(1 + y) 



COS 0i + COS 



2 /i/ 1 



«7 



1 + 1/ 



{l - (1 - y)cos6»icos6»2} (cos 6^ + cos fl 2 ) . (A.19) 



The C B distribution is given by 
dN 



dC B 



oc 



dn*dnidn' 2 \m\ 2 5(c b - c B [e*,^*,^]) 

1 fdAps (1 - y) 2 dN vl y(l - y) dN W2 , «7 d^cvl 



As 



dC B 



l + y dC B 2(1 + y) dC B + l + y dC B / 



(A.20) 



but the last term integrates out to zero: 
dC B 



jj dcos^dcos# 2 (cosfli + cosfla) {1- {l-y)cGsd' 1 cosQ' 2 \ 5(C B - C B [r][,r]' 2 ]) 

/ / di][ dr]' 2 sech 3 rj[ sech 3 r\ 2 sinh (rj[ + rj 2 ) { 1 — (1 — y) tanh r\ x tanh r( 2 } x 
1/ J —00 



x 5 [ C B — tanh 
= 



(A.21) 

Thus, in the VVF case, the chiralities a and 7 do not have an effect on the C B distribution 
at threshold. The second term in ( A.20| ) is 

dAVi 



dC B 



oc J J d cos 6[ d cos 6' 2 cos 2 9[ cos 2 6' 2 5 ( C B — C B [r][ , r/ 2 ^ ) 

/ / 2 / 2 / 2/ 2// 

d^i d?7 2 sech r/ 1 sech r/ 2 tanh 7^ tanh r/ 2 <5 I C B — tanh 



dr]' 2 I cosh — Tj' 2 ^j + l| sech 2 sech 2 r/ 2 tanh 2 r/^ ' 1 tanh 2 r( 2 



l-C B J_ 



sech 2 rj 2 tanh 2 rj 2 tanh 2 ^rj' 2 + In 


"l+C B " 
1-C b 


) 


cosh 2 [r] 2 + In 


r i+c B " 

_1-C B _ 


) 





(A.22) 
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After integrating out rj 2 , we get ( |6.13| ) 



Jvi 



diVyi _ 1 - C\ 
dC B 



24 Cl 



2C B (15 + 8C§ + 10C B + 8C B + 15C B ) 



+ 3 (5 + C\ + 2C B + 2C B + C B + 5C B °) In 



i + e B 
i- c B 



(A.23) 



By a similar calculation, we obtain ( jug ) 



JV2 



d/Vv2 



^ |-2C B (3 + ei - C B - 3C B ) + 3 (1 - Cl) In (I±||) } . (A.24) 



In total, we get the answer in (|6.8|) 

dY 1 ( diVps (1 - y) 2 diVyi y(l - y) dY V2 

l + y dC B " 2(1 + y) dC B 



dC B Y tot 1 dC B 
9(1 +i0 



2(2 + y) 2 1 4(1+ y) 

A. 5 C B distribution in SUSY and MUED 



T (^^ T y(l-S/) T 
2 1 J PS _ , T, , 7 J V1 — . T J V2 



+ 



(A.25) 



For SUSY (MSSM), the threshold behavior of C B is given by ( 6.1C| ), since B is a scalar 
particle (slepton): 



dY 
dC^ 



Jps • 



For the minimal UED model (MUED), the matrix element is 

|A4| 2 oc 1-Ai(cos0i -cos6» 2 ) - A 2 cos 6^ cos# 2 



(A.26) 



(A.27) 



with 



l-2y 
l + 2y 

l-2y 
l + 2y 



o _ M z r 12 



1+3 



(A.28) 



is the mass of the Z boson and Cw = cosOw is the cosine of the Weinberg angle 9w- 
Comparing (A.27) to (A. 5), we find C = —A 2 and from (A. 11) we get 



dY 
dCe" 



jps — t ( - , y \ jf • 



4 \l + 2y 



(A.29) 
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